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• . Abstract 

m 

C^ ■ In this paper, we consider coupled forward-backward stochastic differential equations (FBSDEs in short) 

O^ ' with parameter e > 0, of the following type 

X^-*'- (s) = x + /; / (r, X^-''^ (r) , Y^'*'- (r)) At + ^ j; a {r, X^'''^ (r) , F^'*'- (r)) dW^ (r) , 
Y'^'*'"(s) = /i(A'='*'"(T))-f/Jp(r,X='*'"(r-),F='''"(r-),Z'^'*'"(r))dr 
/\ ■ [ -/JZ'''*'^(r)dH/(r), 0< t< s<T. 

. - -' We study the asymptotic behavior of its solutions and establish a large deviation principle for the corresponding 

processes. 

AMS subject classifications. 60F10; 60H10. 

Key words: Large deviation principle, Forward-backward stochastic differential equations. Random perturbation, 

Meyer- Zheng topology. 

1 Introduction 

Non-linear backward stochastic differential equations (BSDEs in short) were first introduced in Stochastic Optimal 
Control Theory with the pioneering work of Bismut [6] and then developed by Pardoux and Peng [25]. Since then, 
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they have become a powerful tool in many fields, such as mathematics finance, optimal control, stochastic games, 
partial differential equations and homogenization etc. Simultaneously, it is well known that the Hamiltonian system 
associated with studying maximum principle of stochastic optimal control problems is the corresponding to certain 
fully coupled forward backward stochastic differential equations (FBSDEs in short) (see [30]). In mathematical 
finance, fully coupled FBSDEs can be encountered when one studied the problems of hedging options involved 
in a large investor in financial market (see [10, 16[). Besides, fully coupled FBSDEs can provide probabilistic 
interpretations for the solutions of a class of quasilinear parabolic and elliptic PDEs (see [27, 28, 29, 30]). 

There are three main approaches to solve FBSDEs. The first one, the Method of Contraction Mapping, which, 
was considered by Antonelli [2] and later developed by Pardoux-Tang [29], works well when the duration T is 
relatively small. Second, the Four Step Scheme. It removed the restriction on the time duration for Markovian 
FBSDEs, as initiated by Ma-Protter-Yong [23]. Third, the Method of Continuation. This method can treat non- 
Markovian FBSDEs with arbitrary duration, initiated by Hu-Peng ]19] and Peng-Wu ]30], and later developed by 
Ma at al in ]22]. It is worth noting that these three methods do not cover each other. 

In our paper, on one hand, we obtain two kind of asymptotic results. The first employs the assumptions 
established by Peng et al in [19[, [30[ and provides the convergence in distribution of the associated processes 
(Theorem 15). The second, more closely related to the PDF point of view but still using strictly probabilistic 
methods, provides almost-sure convergence with values in some L^ type spaces (Theorem 16). In this last result 
the problem of convergence of classical/viscosity solutions of the quasilinear parabolic system of PDE's associated to 
the backward equation is naturally addressed. Notice that when this PDF takes the form of the backward Burgers 
equation, the problem becomes the convergence of solutions for vanishing viscosity hydrodynamical parameter. 
We mention the work [18[ where the authors do the asymptotic studies of FBSDEs with generalized Burgers type 
nonlinearities. 

The large deviation principle (LDP) characterizes the limiting behavior of probability measure in applied prob- 
ability and largely used in rare events simulation. Recently, there has been a growing literature on studying the 
applications of LDP in finance (see [31[). 

We now consider the following small perturbation of FBSDEs (1.1) 

X^^t^^ {s) = x + ///(r,X^^*'="(r),r'^'*^^(r))dr 

+Ve // cr (r, X^^*'=" (r) , y^'*'^ (r)) dT4^ (r) , 
r^^*'="(s) ^ /t(X'^'*^^(r))+/J.g(r,X^^*'="(r),r^'*'^(r),Z^^*'="(r))dr ^^''^^ 

_ jj z^,^,x (j.) dVF (r) , < i < s < T. 

The solution of this equation is denoted by 

I^X^.t,x ^^^ ^ ye,*,:. (^^ ^ ^cUx (s),t<S<T). 

On the other hand, we want to establish the large deviation principle of the law of y^^*'^ in the space of 
C ([0, T] ; R") , namely the asymptotic estimates of probabilities P {Y^^*''-" e F) , where T e B {C ([0, T] ; R")) . 

Ma et al in [24[ first considered the sample path large deviation principle for the adapted solutions to the 
FBSDEs in the case where the drift term contains Z and under appropriate conditions in terms of a certain type of 
convergence of solutions for the associated quasilinear PDEs. While, with probability methods, by the contraction 
principle, the same small random perturbation for BSDEs and the Freidlin-Wentzell's large deviation estimates in 
C ([0, T] ; R") are also obtained by [13, 17, 18, 32] (see references therein). In ]11], a large deviation principle of the 
Freidlin and Wentzell type under the nonlinear probability for diffusion processes with a small diffusion coefficient 
was obtained. 

Our aim in this paper is extend the previous work ]13, 17, 18, 32] to the coupled case, i.e., when both drift 
and diffusion terms contain Y, by probability methods and under some suitable assumptions. Note that in ]11], 
the authors considered the fully coupled FBSDEs via the corresponding PDF techniques. 



We present the framework as follows. In Section 2, we give the framework of our paper. Some estimates and 
regularity results are established for the solutions of FBSDEs (1.1) in Section 3. Then in Section 4 we show our 
main results Theorem 15 and Theorem 16. Section 5 is devoted to establishing the large deviation results for (1.1). 
Some proofs of technique lemmas are given in the Appendix. 

2 Preliminaries 

Let us begin by introducing the setting for the stochastic differential systems we want to investigate. Consider as 
Brownian motion W is the d-dimensional coordinate process on the classical Wiener space {ft,T, P), i.e., ft is the 
set of continuous functions from [0, T] to R'' starting from 0, Q ~ C {\0,T]; R'^, O) , J^ the completed Borel cr-algebra 
over fi, P the Wiener measure and W the canonical process: Ws (oj) = oJs, s G [0,T] , to ^ ft. By {J^J, < s < T} 
we denote the natural filtration generated by {M^s}o<s<t ^^'^ augmented by all P-nuU sets, i.e., 

T,^a{W{r),r<s}\/Mp, se [0,r], 

where Afp is the set of all P-null subsets. For each i > 0, we denote by {J^*, t < s < T} the natural filtration of 
the Brownian motion {W (s) — W (t) , t < s < T} , augmented by the P-null set of T. 
We denote by Ai^ (t, T; Rp) the set of all RP-valued J't-adapted process v (•) such that 



E 



\v {s)\ ds 



< +00, 



and S^ (t, T; R") , the set of all M"-valued J't-adapted process v (•) such that 

E 



sup \v {s)\ 

t<s<T 



< +00. 



Clearly, TVt = S^ {t, T; M") x S^ (t, T; R") x M^ (t, T; Rp) forms a Banach space. 

The space Aft is a Banach space for the natural norm associated with the product topology structure. 
We use the usual inner product and Euclidean norm in R", R™, and M™^''. 

Definition 1. A triple of processes 

(X^'*'^,r^'*'^,Z'^'*'^) : [0,T] X SI ^ R" X M" X R"^'' 

is called an adapted solution of the Eqs. (1.1), i/ (X^^*'^, F^'*'^, Z^'*'^) G A/*, and it satisfies (1.1) P-a.s.. 

It is clear that the above Eq. (1.1) is a stochastic two point boundary value problem. Especially, it contains 
deterministic two point boundary value problem as a special case when e — > 0. 
We are given an tti x n full-rank matrix G. We use the notations 

A''{t,u) = 





We now give the first assumptions of our paper: 

(Al) Let m > n. The processes g {-jX^y, z) , / (•, x, y) and cr (•, a;, y) are J^t-adapted, and the random variable h (x) 
is J^T-nieasurable, for all {x,y,z) G R" x R™ x R™^*^. /, a and g are continuous in t, P-a.s.. Moreover, the 



following holds: 

|2' 



E 



\h{OW 



< +00. 



There exists a constant Ci > 0, such that, u* = {x\y\ z') e K" x K™ x M"^'', i = 1, 2, 

\g{t,x^,y^,z'^) - g(t,x^,y^,z^)\ < Ci {\x^ - x^\ + \y^ - y^\ + \z^ - z^\) , 
\f{t,x\y^) - f {t,x^,y')\ < C\ (|xi -x2| + \y^ - y^\) , 
\a {t, x\y^) - a {t, x^,y^) \ < Ci {\x^ - x^\ + \y^ - y^\) , 
\h{x^) -h{x^)\ <Ci\x^-x^\, 
P-a.s., a.e. t e M+. 

(A2) There exists a constant C2 > 0, such that for any fixed e E (0, 1] , 

{A'{t,u')-A'{t,u^),u'-u^) 
< - (C2 + Vi) (\G {x^ - x^) f + |G^ (yi - y^)\^) , P-a.s., a.e. t e M+, 

and 

{h{x^) ^h{x^),x^ -x^) >C2\G{x'^ -x^)f , P-a.s., V(a;\x2) e M" x R". 

We have the following: 

Proposition 2. Assume that (Al) and (A2) hold, then there exists a unique adapted solution (x^'*'^,y^'*'^, Z^'*'^) 
Jor Eqs. (1.1). 

The proof can be seen in Theorem 2.6, Remark 2.8 in [30]. 

Remark 3. If m = n, G ^ /„. If m < n, according to Theorem 2.6 in [30], the first condition of (A2) should be 

{A'{t,u^)~A'{ty)y~u'') 

However, if setting x^ ~ x^, y^ ~ j/^, there exists a contradiction, because a does not contain z. 
From now on, let us suppose that 

the coefficients /, a, g and h are deterministic, i.e., independent of cj G J7. 

Now consider 

m'^ (i, x) = r^'*'^ (i) , (t, x) € [0, T] X M, e G (0, 1] , (2.1) 

which is a deterministic vector since it is J^f measurable ( Blumenthal's 0-1 Law -see remark 1.2 of [14]). In ]28] it 
is shown that u is a viscosity solution of the associated quasilinear parabolic partial differential equation 

-Q^—'y^^ ^>+2 ^*=i ""^ ^ ' ^' ^" ' ^'^ dxdx- ' ^' ^*=i ■''^ ' ^' ^^ ' ^^' dx ^ ' ^' 

+g\t, X, u^it, x), ^V^u^{t, x)a^t, x, u^t, x))) = 0, ^'^''^' 

y{T,x) = h{x), xG M"; tG [0,r]; l^ l,...,n, 

where aij — {acr'^)ij. 



We define below the notion of viscosity solution for the parabolic system of Partial Differential Equations (PDEs 
for short) (2.2). For each e > 0, consider the following differential operator, 

/ = l,...,n, V^eC^'2([0,r] xM",M"), te [0,r], (x,2/,z)e]R" x K" xM"^''. 

The space C^-'^{[0,T] x R",]R") is the space of the functions (p : [0,T] x R" ^ M", which are C^ with respect to 
the first variable and C^ with respect to the second variable. 
The system (2.2) reads 

5' {t, X, u^t, x), ^/eV^u''{t, x)a{t, x, u" (t, x))) = 0, (2-^) 

w^(T, x) = h{x), X e M'', t e [0, T], / = 1, ..., n. 

Definition 4 (Viscosity solution). Let u^ G C([0,r] x M",]R"). The function u^ is said to be a viscosity 
sub-solution (resp. super- solution) of the system (2.2) if 

{u'y{T,x)<h\x); yi^l,...,n; a; e M" 

(resp. (a^)\T,x)>h\x); V/ = l,...,n; x <eW ) 

and for each I — 1, ...,n, (t,x) G [0,T] x M", (^ G C"'"'^([0,r] x ]R",]R") smc/i that {t,x) is a local minimum (resp. 
maximum) point of Lp — u^ , we have 

-^{t, x) + {Llip)(t, X, u^'^t, x), ^\/a;ip{t, x)a{t, X, u" (i, x))) + 
g^(t,x,u'^{t,x), y/eWxf{t, x)a{t,x, u^(t,x))) > 0;l ~ 1, ...,n 

(resp. < 0). 

The function u^ is said to be a viscosity solution of the system (2.2) if u^ is both a viscosity sub-solution and a 
viscosity super- solution of this system,. 

Under more restrictive assumptions (that we present in Section 3), we shall have 

us(i,x)=y*'"'^(i) (2.4) 

and w'^ will be actually a classical solution of (2.2). This can be proved using the Four Step Scheme Methodology 
of Ma-Protter- Young [23] with the help of Ladyzhenskaja 's work in quasilinear parabolic PDEs [20] . 

For the simplicity of notations, we only consider the case where both X^^*'^ and y^^*'=^ are n-dimensional, that 
is TO = n, G = /„. The result is analogous for the case m > n. 

Next we introduce another set of assumptions which is slightly different from (A1)-(A2). 

(A3) We say /, g, a, and h satisfy (A3) if there exists two constants Ci, A > such that : (Al) hold as well as: 

(A.3.1) yt e [0,T], V(a;i,yi,zi), (x2,y2,Z2) G M" x M" x K"^^: 

{xi-X2,f{t,xi,yi)-f{t,X2,yi)) < Ci\xi-X2\ , 
{yi-y2,f{t,xi,yi)- f(t,x,iy2)) < Ci\yi-y2\^, 



(A.3.2)yte [0,T],y{x,y,z) e K" x M" x 



pnxd. 



\f{t,x,y)\ < A(l + |x| + |y|), 

\git,x,y,z)\ < A(l + |a;| + |y| + |z|), 

\a{t,x,y)\ < A(l + |a;| + |y|), 

IM^)I < A(l + N). 



(A.3.3)yte [0,T],y{x,y,z) e K" x M" x 



ixd. 



u H^ f{t,u,y), 

V H- > (7(t, X, w, z) are continuous mappings. 

Under this set of hypothesis, theorem 1.1 of [14] ensures that there exists a constant C — C{L) > , depending 
only on Ci, such that for every T < C (1-1) admits a unique solution in Mt- 
Moreover, using theorem 2.6 of [14] , we have 

P(Vs e [t,T] : u%s,X'/)^Y',-^) = 1. (2.5) 

P®^({(u;,s)ef^x [i,T] : |Z*''^'^(u;)| > Ti}) =0. (2.6) 

where n stands for the Lebesgue measure in the real line and Fi is a constant which only depends on Ci, A, n, d, T. 

By remark 2.7 of ]14] we have that, for each e > 0, u*^ only depends on the coefficients of the system (2.1), 
f,g,^/ea,h. The fact that the dependence of Fi — Fi (Ci, A, n,T) determines that the properties (2.5), (2.6) 
above hold uniformly in e; in particular, there exist continuous versions of (Y^'^'^ , Zl'^''^)t<s<T which are uniformly 
bounded, also in e. 

We now assume more regularity on the coefficients of (1.1) in the next set of assumptions. 

(A4) For T < C, we say that /, 5, /i, a satisfy (A4) if there exists three constants A, A 7 > 

(A4.I) yte [0,T],y{x,y,z) G R" X M" X R"^'': 

\f{t,x,y)\ < A{l + \y\), 

\git,x,y,z)\ < A(l + |y| + |z|), 

\a{t,x,y)\ < A, 

\h{x)\ < A, 

(A.l2)y{t,x,y) e [0,r] X M" X R" X R"^^'' : {^,a{t,x,y)^) > A |C|\ V^ G R", where a(i,x,2/) ^ aa'^ {t,x,y). 

(A. 4-3) the function a is continuous. 

(A. 4-4) The function a is differentiable with respect to x and y and its derivatives with respect to x and y 
are 7- Holder in x and y, uniformly in t. 

Under the set of assumptions (A. 3) and (A. 4), using proposition 2.4 and proposition B.6 of [14] and theorem 
2.9 of [15], one can prove that there exist two constants, k,ki, only depending on i. A, n,r (independent of e) 
such that : 

\u'{t,x)\<K, (2.7) 

u' eCl'^{[0,T]xR"), (2.8) 



sup \VxU%t,x)\< Ki (2.9) 

{t,x)e[0,T]xR'' 



and 



^t^^e.x ^ ^V^u'^(i,X*'^'^)cr(s,X*^%y,*^^'=^). (2.10) 

where u^ solves uniquely (2.2) in C^' ([0,T] x R"), the space of C^ functions with respect to the first variable and 
C^ with respect to the second variable, with bounded derivatives. All these facts can be proven probabilistically. 
The last claim and the properties (2.7), (2.8) and (2.10) are proved in [14]. Delarue delivers in the appendices of 
[14] probabilistic methods to obtain these regularity results, under assumptions (A. 3) and (A. 4) and over a small 
time enough duration, using Malliavin Calculus techniques. The estimate of the gradient (2.9) is established by 
a probabilistic scheme in [15], using a variant of the Malliavin- Bismut integration by parts formula proposed by 
Thalmaier [33] and applied in [34] to establish a gradient estimate of interior type for the solutions of a linear 
elliptic equation on a manifold. 

To start with, let us fix x, y e R" , £ G (Oj 1] • We establish second order moment estimates for the solution of 
FBSDEs (1.1), which will be essential in Section 3. For convenience, we use the following notations in this paper: 



f /^'*'-(r) 
(r) 



a 



gS,t,X J^j,-| 



/(r,X='*'^(r),y^'*'^(r)), 

cr(r,A:^'*'^(r),r^'*'^(r)), 

g (r, A:^'*'^ (r) , F^'*'^ (r) , Z^'*'^ (r)) 

/i (A-^'*^^ (T)) , 



etc. 



The follwing proof will be found in Appendix 

Lemma 5. Assume that (Al) and (A2) hold. Then we have 



E 


sup |A^'*'^(s)-A'^'*'!'(s)|^ 

t<s<T 


< Ci\x-y\\ 


E 


sup |y^'*'^(s)-y"'*'2'(s)|^ 

t<s<T 


< Ci\x-y\\ 


E 


/f |Z'^'*'^(s)-Z^'*'2'(s)|^rfs 


< C,\x-y\\ 



where Ci is a positive constant independent of e and t. 
Lemma 6. Assume that (Al) and (A2) hold. Then we have 



E 


sup |A^'*'^ (s)r 

t<s<T 


< 


C2{1 


E 


sup |r'^'*^^(s)|^ 

t<s<T 


< 


C2(l 


E 


'j^\Z^-^^-^{sfds 


< 


C2(l 



where C2 is a positive constant independent of e and t. 
The following Lemma shows the continuity on t. 



(2.11) 



(2.12) 



< C3\h-t2\{l + \x 



< Calti-tal l + lx 



Lemma 7. Assume that (Al) and (A2) hold, then we have. 

E sup |X^'*i'^(s)-X'^'*2^^(s)| 

.tlVt2<S<T 

sup |y^'*i'==(s)-y^'*2'="(s)r 

.tiVi2<s<T 

where C4 is a positive constant independent of e. 

Lemma 8. Assume that (Al) and (A2) hold. Pick < £2 < £1 < 1. Then we have 



< Calti-tsl 1 + b 



E 



sup |X^i'*'^(s)-X^2'*'^(s)r 

t<s<T 

sup |r=i'*'^(s)-y'^2'*'^(s)|^ 

t<s<T 

J^ iZ^ut,'^ (s) - Z^^'*-'' {s)f ds 



< C4(veT-V^) 



where C3 is a positive constant independent of e and t. 
Now consider the following deterministic equations 

J A'*'^(s) = a; + // / (r, A-*'^ (r) , y '^ (r) , 0) dr, 



(2.13) 



(2.14) 



(2.15) 



[ y'^(s) = /i(A:'*'^(r))+// g(r,A'*'^(r),y'^(r),0)dr. 

By the regularity of /, g, h we have the following 

Lemma 9. Assume that (Al) and (A2) hold, then there exists a unique solution (A"*'^, 3^*'^) for Eqs. (2.15). 
Clearly, the estimates in Lemma 8 shows that the triple of 

(x^'*'-(5),y^'*'-(5),z^-*'-(s))^^j^^^, 

are Cauchy sequences and therefore converge in 

(52 (i, T; R) X 52 {t, T; M) X A^2 (^^^ y . j^d^ ^ ^ 

Denote the limit by (A'"'*'^ (s) , 3^0,*, a; (^g^ ^ q^ . , as e ^ 0. By uniqueness and existence of Eqs. (2.15), we know 

that the limit (A''''*'^ (•) ,3^0'*'^ (•) ,0) is the unique solution of Eqs. (2.15). Therefore, the conclusions in Lemma 
5, Lemma 6, and Lemma 7 also hold when e — > 0. 



3 Main Results 

3.1 Convergence of distributions 

In this subsection, we first study FBSDEs (1.1) with small noise intensity. To conclude this subsection we introduce 
the notions of pseudo-path topology and quasimartingales (Meyer- Zheng [21]), adjusted to our setting. For that, 
let us introduce the following notations. 

Let T > be a real constant (the terminal time). For any natural number I > Q,hy D (K.') = D ([0,T] ;]R') 
we denote the Skorohod space of right continuous with left-hand limits functions a; on [0, T] with values in M' such 



that 



X (T— ) = limx (t) = X (T) 



and by convention, x (0—) = 0. Additionally, we introduce the following metric p on D (R') : 

p{x,y)^ [ {\x{.s)-y{.s)\Al)d.s, x,yeD{R'). 



The topology induced by this metric is the Meyer-Zheng topology introduced below on D (R') . Let ( be the 
coordinate mapping on D (R') defined by 

Ct{x)^x{t), te[0,T], xeDiR""), 

and introduce the c-algebras of subsets of D (R') , 

Vl ^ Vl (R') =a{{Cu:ue [t, s]}) , < t < s < T, 

V = V (m') = 2??, (r') . 

(Meyer-Zheng topology) Let A (dt) the measure e^*dt on M+. Let w (t) be a real valued Borel function 
on R_|-. The pseudo-path of w is a probability law on [0,oo] x R : the image measure of A under the mapping 
t — > (t, w(t)) . We denote by ip the mapping which associates to a path w its pseudo-path: it is clear that ip 
identifies two paths if and only if they are equal a.e. in Lebesgue's sense. In particular, V' is 1 — 1 on D (R') , and 

provides an imbedding of D (R') into the compact space V of all probability laws on the compact space [0, oo] x R . 
We give to the induced topology on D (R') the name of pseudo-path topology or Meyer-Zheng topology. Let us 
introduce some intermediate sets between D (R') and V: ^ will be the set of all pseudo-path. We have inclusions 

D (r') C * C ^. 

The following characterization of the Meyer-Zheng topology is worth noting. 

Lemma 10. The pseudo-path topology on ^ is equivalent to the convergence in measure. 

Furthermore, it is known that ^ is a Polish space; and D (R') is a Borel set in V. 

Lemma 11. Let B (D (M!')) be the a-algehra of Borel subsets of D (M!') in the Meyer-Zheng topology. Then 
B{D{M.^)) =X>(R'). 

The proof can seen in [21]. The most important application of the Meyer-Zheng topology is a tightness result 
for quasimartingales. We give the definition here 

Definition 12. Let X be an J- -adapted, cadlag process defined on [0, T] , such that E\X [t)\ < oo for all t > 0. For 
any partition tt : = tg < ti < • • • < t„ < T, let us define 

Vf{X):= Y, E{|IE{X(t,+i)-X(t,)|-^tJ|} + E|X(i„)|, (3.1.1) 

0<i<n 

and define the conditional variation of X by Vt (X) := supV^ (X) . If Vt {X) is finite, then X is called a quasi- 

TT 

martingale. 

The following result holds. 



Proposition 13. Let P„ be a sequence of probability laws on D ("R') , such that under each Pn the coordinate 
process X {■) is a quasiniartingale with conditional variation Vn {X (•)) uniformly bounded in n. Then there exists 
a subsequence (Pn^) which converges weakly on D (M') to a law P, and X {■]}] is a quasiniartingale under P. 

The proof can be seen in [21]. 
Now we turn back to Eqs. (1.1). Define 



u^{t,x) = Y^-*^^{t), 
v^ (t, x) = Z^'*'^ (t) , 



etxn^ for {t,x)e [0,T]xM. (3.1.2) 



From the existence and uniqueness of solution for Eq. (1.1), we have the following Markov property 

\v^{s,X^'^--{s)) = Z--^^^^'''^(^) (5) = Z^^*--(5) a-e.a.e. se[t,JJ. (3.1.3) 

In Lemma 5, we have already seen that the function u^ is Lipschitz continuous in x, uniformly in t. With the help 
of Lemma 7 we also know the continuity property for u'^ in t. We now set 

fis^x) = f{s,x,u''{s,x)), 

a'^ {s,x) — a {s,x,u^ {s,x)) , (3.1.4) 

g^s,x) = g{s,x,uUs,x),v'' {s,x)). 

We study the properties of /^, ct^, p^ and ct^ as follows: 

Lemma 14. Assume that (Al) and (A2) hold. Then f^, a^ /°, and a^ satisfy uniformly Lipschitz continuous and 
/^, ij^ converge uniformly /°, a'^ , respectively. Moreover, f^, a^ , f^ and a'^ satisfy sublinear growth. 

Proof For any x, y G M", s G [i, T] , we have 

If {s,x)- f {s,y)\ = \f {s,x,u'' {s,x))- f {s,y,u^s,y))\ 

< Ci{\x-y\ + \u^s,x)-u^s,y)\) 
= Ci\x-y\ + CiCi\x-y\ 

< max{Ci,CiCi} |a; - y| , 

by Lemma 5 and (Al). The same properties for a^ f^, and /" are proved similarly. Next, we show the uniform 
convergence, 

\r {s,x)- f^{s,x)\ = \f {s.x.u'' {s,x))- f [s,x,u° {s,x))\ 

< Ci\u''{s,x)-u°{s,x)\<CiC3V£, 

by Lemma 8 and (Al). Once again, the same properties for a^ are proved similarly. We are going to prove b"^ 



^Note that the quasiniartingale in [21] is defined on [0, +oo]. However, it is fairly easy to check that if X is a quasiniartingale on 
[0, T] as is defined above, then the process X (t) = X (t) I[o,T) (*) + ^ (^) I[T,+oo) (*) , * S [0, +oo] is a quasiniartingale in the sense 
of [21], Furthermore, the conditional variation Vt (X) defined here, is exactly the same as V {X) defined in [21], In other words, our 
quasiniartingale is a "local" version of the one in [21]. 
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satisfies sublinear growth. We have 



|r(s,x)-/(s, 0,0)1 = \fis,x,u'is,x))-fis,0,0) 

< C,{\x\ + \u'{s,x)\) 

< Ci(|x| + VC^(l + |x|)) 



The assumption (Al) yields 



\r{s,x)\<(c\+C\^)\x\+C\^+ sup |/(r,0,0)| 



0<r<T 

Once again, the same properties for a^ , /", a'^ are proved similarly. 

Then it is easy to check that {X"-*'^ {■) ,Y^-*^'' (•) ,2'=^*'^ (•)) solves the following decoupled FBSDEs 

X^^t^^(s) = a; + ///'^(r,X^^*'^(r))dr + V£//cr^(r,X'^'*'^(r))dM^(r), 
r^.*^^ (s) = h (X"^'*'^ (T)) + /J g^ (r, X^^*'"' (r)) dr - /J Z^'*^^ (r) dVK (r) , 
0<t<s<T. 

From now on, we are concern on the phenomenon about the laws of (X^'*'^, y^'*'^) when e — > 0. 

Theorem 15. Under the assumptions (Al) and (A2), we can conclude the following results: 
i) For all (5 > 0, 

lim P I sup |X"^*'=" (s) - X^'"" {s)\ xjI = 0. 

e^" \t<s<T J 



D 



(3.1.5) 



(3.1.6) 



ii) Let Q'^ ^ P I (y^^*'^ (•)) ) he the probability measure on D(]R"). Then there exists a subsequences Q^" 
of Q^ and a probability law Q on D (K") such that Q"^" converges weakly in the Meyer-Zheng topology to Q as 



n — )• +00. 



Proof. It follows from the definition of X^' '^ and X '^ that 

sup |X'^'*^^(s)-A'*^^(s)| < / |/"(r,X^'*^^(r))-/"(r,A'*^^(r))|dr 

t<s<T Jt 



-\fe sup 

t<s<T 



fT(r,X^'*'^(r))dPF(r) 



(3.1.7) 
From Chebyshev's inequality and the first assertion of the Lemma 14 we obtained an estimate of the first term of 
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the right side of (3.1.7): 

P{Jl \r {r,X^'*^- (r)) - f {r,X'-- (r))| dr > |} 
I// \r {r, X^'*-- (r)) - f (r, X'-- (r))\ dr\ 



<4J-2 

= 46-^TE 



< 



^TE 



// \f^ (r, X^'*'^ (r)) - f {r, A"*'^ (r)) | dr 

// If (r, X^'*'^ (r)) - f" (r, <Y*'^ (r)) + /^ (r, A"*'^ (r)) - f (r, A"*'^ (r))|^ dr 

/f I/- (r,X--*'- (r)) - /- (r, A-*^- (r))|VTC3e 



(3.1.8) 



<85-2y(max{Ci,CiCi})^E // IX-^'*'^ (r) - A"*'^ (r)r dr + T (CiCg)^ e 
< 8(5-2t (max {Ci, Cid})^ TE [Cs^i + (CiCg)^ e 

The estimation of the second term in (3.1.7) can be accomplished with the use of the generalized Kolmogorov 
inequality for stochastic integrals: 



p\y/e sup |rcr'^(r,X^'*'^(r))dVF(r)| > ll 

L t<s<T J 



< 4S-^sE 
= 4S-^eE 

< iS-^eE 

< 4S-^s 



/; |cr'^(r,X'^'*^^(r))rdr 

/f |cr (r, X^'*'^ (r) , r^^*'=^ (r)) - a (r, 0, 0) + a (r, 0, 
4Ci2 /f (|X^'*'^ (r)l' + |r^'*'^ (r)l' + 2 |(T (r, 0, 0)|') dr 
8rC2C2 (l + \x\') + J^ 2 \a (r, 0, 0)|' dr 



' dr 



(3.1.9) 



by Lemma 6. Estimates (3.1.8)-(3.1.9) imply the assertion of the theorem (i). 

We are going to prove the second one. First, we establish the connection of solution of BSDE (3.1.5) to 
quasimartingales. Given a subdivision tt : = to < ^i < ■ ■ ■ ^n = 2^, we get 



v^ (y^'*'^) 



n-l 

= E\h{x'-'^-{T))\+ EE[|E[r^'*'-(tfe+i)-r^-*'-(tfe)|j-t]|] 

fe=0 

"-If { t 
= E\h (X-^'*'^ (r))| + X; E E / "+' g (r, X^'*'^ (r) , F^'*'^ (r) , Z^'*'^ (r)) dr 

< E\h (X-^'*'^ (r))| + E [/q^ g (r, X^'*'^ (r) , r^'*'"= (r) , Z^'*'^ (r)) dr 

< E|/i(x'^'*'^(r))-;i(o) + /i(o)| 

/f |5"'*'" (0 - 3 {r. 0, 0, 0) + 5 (r, 0, 0, 0)| dr 

< CiElX-^'*'^ (T)| + CiE [/f dX-^'*'^ (r)| + [y^'*'^ (r)| + jZ-^'*'^ (r)|)dr 
+ /f |5(r,0,0,0)|dr + /^(0). 



J"t 



(3.1.10) 



By Jensen's inequality, it follows from Lemma 6 



E|X'^'*'^(T)| < WCa 1 
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and 



ft i\X^'*''' {r)\ + |y^'*^^(r)| + |Z^'*^^(r)|)dr 



(e/^ |X^-*'^(r)|^dr)' + (e/^ |r^'*^^(r)|^dr)' + (e/^ |Z^'*'^ (r)|^ dr 



< \/TE 

< Vt 
<3Vt{tc2{i 

by Holder inequality and Jensen's inequality for concave function. So we have 



y. (ye.t,x) < ^ /^^ /^ ^ |^|2\ ^ 3^ /^^ A 



+ /i(0)+T sup |5(r,0,0,0)|. 

0<r<T 

Hence, by (Al), noting that Vt (X) {Y^'*^^) = supK^ (ye,*,^) < +00, the resuh follows. 

Now since D (M") is a separable metric space, there exists a compact metric space K such that D (M") is a 
subset of iiT. Note that D (K") is a Lusin space: for every embedding in a compact metric space K, D (K") is a 
Borel set in K: D (R") € B (K). On the compact metric space K we define 

Q^ (^) = 0" (A n Z) (M")) , ^e6(/f). 

Clearly, Ar\D{W) belongs to S {D (K")) = X> (M"), the last equality being true in view of Lemma 11. The set of 
probability measure on the compact metric space K is compact in the weak convergence. Hence, we can choose a 
subsequence also denoted (£n)„>i and a probability measure Q on K such that 



We now show that 



Q"" ^QonK. 



Q(i5(M")) = l. 



We notice that Q^ is the distribution of Y^'^ considered as a random variable with values in {K, B {K)) . Further- 
more, by Proposition 13, we know that, possibly along a subsequence, Q^" converges weakly to a probability law 
Q* G Ai {D (K.")) . The uniqueness of the weak limit implies that 



In particular, 

The proof is complete. 



l = g*(i?(M")) = Q(i?(K")). 



(3.1.11) 



D 



3.2 Almost sure convergence 

Considering assumptions (A3)-(A4), we have the following 

Theorem 16. Under the assumptions (A3), and (A4), assuming further T < C, where C = C'{Ci,T) depending 

on the Lipschitz constant Ci and on T , we have 

1. For each s,t £ [0, T], t < s, the solution of (2.1), [X^'^'^ (s) , F^'*-^ (s) , ^e^*-^ {s))t<s<T converges in Aft, when 
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e — > to {X (s) , y (s) , 0)t<s<T7 where {X (s) , Y {s))t<s<T solves the coupled system of differential equations: 

(3.2.1) 



Xis)=fis,Xis),Yis)), 

Yis)^ = -gis,X{s),Y{s),0), t<s<T, 

X{t) ^x, Y{T) = h{X{T)). 



2. Denoting u{t,x) = y*'^ {t) , the limit in e ^>- ofY"^'*'^ (t), the function u is a viscosity solution of 

-Q^+Et=iMt,^,u{t,x)) — it,x)+g'{t,x,u{t,x),0)^0, 
y{t,x) = /i(x); x<EW^]t<E [0,T]; / = l,...,n. 



(3.2.2) 



3. The function u is hounded, continuous Lipschitz in x and uniformly continuous in time. 

4- Furthermore, ifu^Cf^' ([0,T] xK"), since (3.2.1) has a unique continuous solution, the function u is a classical 

solution of (3.2.2). 

Proof. Given e, ei > 0, x G M", for T < C, if (X|^*'=", r/^*'=", Z*^^)t<,<T is the unique solution in Mt of 

' j^M,x == X + // /(r, X^^*'=" (r) , r-^'*^^ (r))dr + ^e // cr(r, X-^'*^^ (r) , y^^*'=" (r))dVK(r), 
y/'*'^ == hiX"'^-'-^ (T)) + /J g(r, X"^'*'^ (r) , y^'*'^ (r) , Z^'*'^ (r))dr - /J Z-^'*-^ (r) dVF(r), (3.2.3) 

a; e K''; < t < s < T, 

and (Xfi'*'^, y/ 1'*'"^, Z^^-'^'^)t<s<T is the unique solution in A/t of 

'Xf '*'^ = .T + // /(r,X'^i'*'^ (r) ,y^i'*'^)dr + Vel// cr(r,X=i'*'^ (-^) ^yei,t,x (^r))dWir), 

a; e R'^; < i < s < T. 

Write, for sake of simplicity on the notations, with 5 — e, ei 

f\s) = f{s,X'^''-{s),Y'''^^{s)), 

g'is) = 5(5,X^'*'-(s),y^-*'-(s),Z^'*'-(s)), 

g\s) = ^^5a(s,X^''^^{s),Y^''^^{s)), 

h^ = hiX^^*'"" (T)). 

As in the proof of theorem 1.2 of [14], using Ito's Formula: 



E 



sup \X' {s) ^ X'' {s)Y 

,t<s<T 



< E 



+2E 
+2E 



\a%s)-a''{s)\^ds 



sup / {X'{r)-X'^{r),r{r)-f'^{r))dr 

t<s<TJt 



sup / {X''{r)-X^^r),{a%r)-a^'{r))dW{r)) 

t<s<TJt 



(3.2.5) 
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Using Burkholder- Davis- Gundy 's inequalities^ there exists 7 > such that: 



E 



sup |X^ (s) - X"^ (s) 

t<s<T 



< E 



+2E 



-27E 



|cr'^(s)-cr^i(s)rds 



sup / {X^ [r) - X^^ (r) , /^ (r) - /^^ (r)) dr 

tKsKTJt 

1/2 

|X'^ (r) - X'' (r)l^ |cr"(r) - cr-^HOI^ dr 



(3.2.6) 

where 7 > only depends on the Lipschitz constant Ci and T. Using the Lipschitz property (A.l) , there exists 
7 > 0, eventuaUy different, but only dependent on Ci,T such that: 



E 



sup \X' {s) - X'' {s)f 



t<s<T 



< 



j{eJ^ {\X^ is) - X^^ {s)f + \Y^ is) Y^^ {s)\')ds 

+E \a'(s)-a''{s)\^ds + 2\V^-^\^ 
+2Vt sup \X^s)-X''{s)f\. 

t<s<T J 



(3.2.7) 



Then, assuming that 1 - 2jVT > 



(1 - 2jVt)E 



sup \X^ (s) - X^i (s) 

,t<s<T 



< 



j{eJ^ {\X^ (s) - X^^ {s)f + \Y^ (s) Y^^ {s)\') 
+2 I Vi - V^l' + E / \a%s) - a^i (s)|' d.s}. 



(3.2.8) 



Using estimates as before, we get for some new 7 > eventually: 



E 



< 7 



sup \Y^ (s) - r^i (s) 

,t<s<T 



■E 



\Z' (s) - Z'' {s)\^ ds 



E|/i"-/i"i|Ve( sup / {Y^{s)-Y^'{s),g'{s)-g^'is))ds) 



Using (3.2.7) and the assumption (A. 3) modifying 7 eventually: 
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sup \X' (s) - X'^ {s)f + sup \Y'{s)-Y''^{s)f+[ \Z' (s) - Z'^ {s)f ds 

'<s<T t<8<T Jt 

< j{E\h' -h'^^+E f \X'{s)~X'^{s)\^ds + E f \Y'{s)^Y'^{s)fds 

+E I |a^(s)-a^i(s)|'ds+|Vi-V^f 

+E f \Z^- (s) Z^^ is)\ (|y- (s) Y^^ is)\ + \X^ (s) X^^ is)\)ds}. 



So, there exist C* < C and F, only depending on Ci and 71, such that, for T — i < C* 

E snp \X'{s)-X''{s)f+ sup \Y'{s)-Y''{s)f+ \Z' {s) - Z'' {s)\^ ds 

t<s<T t<8<T Jt 

< -f{E\h'-h''f+E \a'{s)~a''{s)fds+\^-^\'' 
+E( / \Y' is) Y'^ (s)| + \f{s) g'^{s)\ ds)'}. 



(3.2.9) 

We can repeat the argument in [T — 2C*;T — C*] and recurrently we get this important a-priori estimate (3.2.9) 
valid for all [i,r], with some (possibly different) constant 7 > only depending on Ci. 
Now, from (A. 3), the following estimate holds: 



\h' - h' 



< L-'E 



sup \X' (s) - X^' (s) 

t<s<T 



(3.2.10) 



One can also derive the estimate 



J (|ViT^(s)-V^a^ns)r)ds 



< ci IVe- \/ei| + C2IE 



{\X^s)-X^^is)\' + \Y^s)-Y^^{s)\')ds 



(3.2.11) 



where ci, C2 > only depending on T,L and on the bound of \a\. 

Furthermore, using Jensen's inequality and the Lipschitz property of g, for some C3 modified eventually along 
the various steps. 
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E 



\YUs)-Y-Us)\ + \g^{s)~g-^{s)\ds 



< E / \Y- is) - Y-^ {s)\' + 1/ is) - f^ {.s)\' + 2 \Y- (s) ~ Y^^ (s)| I5- (s) - f^ (s)| ds 



< C3 

< CgJE 



\X^ (s) X^^ {s)\' + \Y^ is) - Y^^ {s)\' + \Z^ is) - Z^^ {s)\' 



sup \X' (s) - X"' (s) 

t<s<T 



■E 



sup |y^ (s) - Y"^ (s) 

t<s<T 



sup |Z'' (s) - Z"^ (s) 

t<s<T 



Using (2.10) and the boundedness of \a\^ we have E /^ \Z^ (s)| ds 



as e ^- and 



hm I 

[e-£iKO 



|Z^(s)-Z'^i(s)rds 



}■ 

(3.2.12) 
(3.2.13) 



Furthermore, by Burkholder-Davis-Gundy's inequalities, there exists a constant 7 > 0, eventually new, such that: 



sup \X^s) ^ X^' {s)f + sup \Y^s) ~ Y'^ {s)f + \Z^s) - Z^' is)f ds 

t<s<T t<s<T Jt 



< 7E 

< 7E 



\X^ is) - X^^ (s)\' + \Y^ (s) Y^^ (s)\' + \Z^ is) Z^^ {s)\' ds 



sup IX" (r) - X^' (r)r + sup {Y" (r) -- Y"^ (r)r + sup {Z" (r) - Z"' (r)r ds 

t t<.r<.s t<.r<.s t<.r<.s 



(3.2.14) 



Moreover, for some 71,72 > , using (3.2.10)- (3.2.14) 



E 



sup \X' (s) - X^' {s)\^ + sup \Y'{s)-Y''{s)\^+ I \Z^s) - Z^' {s)\^ ds 



t<s<T 



t<s<T 



< 7l|^/^-\/el| +72 1 



/ sup \X' (r) - X^' (r) 

Jt t<r<s 



+ sup |r^(r)-r^i (r)r+ sup \Z^r) - Z^^ {r)\'' ds 

t<r<s t<r<s 



(3.2.15) 
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Using Gronwall's inequality, 



E 



/ sup IX" (r) - X"' (r)l^ + sup {Y" (s) - Y"' {s)\^ dr 

Jt t<r<s t<r<s 



< CalV^-V^I 



/ sup |Z''(s)-Z^i (s)l^ds 

Jt t<s<T 



(3.2.16) 
is bounded, by 



for some C2 > only depending on Ci (independent of e), since E sup(<j,<y \Z'^ [s) — Z^^ (s) 
the results (2.9) and (2.10) . 

We conclude, by the previous estimates, that the pair (X*'^ (s) ,y*'^ (s))t<s<T form a Cauchy sequence and 
therefore converges in 5^ (t, T; W^) x 5^ (t, T; K*^) . 

Denote (X (s) , Y (s))t<,<T its limit. {X (5) , y (s) , 0)t<,<T is the limit in^ of (X*'- (s) , F*'^ (5) , Z*'^'- (s))t<«<T 
when e ^ 0. 

If we consider the forward equation in (1.1) and if we take the limit pointwise when e — > 0, we have 

X*{s)=x+ f f{r,X{r),Y{r))dr, (3.2.17) 

where we have used the boundedness of a and the continuity of /. 

Similarly, we can take the limit on the backward equation when e — > 0. Using the continuity of the functions 
h,g and 



E 



Z*'^ (r) dT4^(r) 



E 



Z'^'^ir)] dr 



as £ — > 0, 



which implies J Z'^'^ (r) dTy(r) -^ 0, a.s P, we have 



n* = h{X (T)) + / g{r, X{r),Y (r) , 0)dr 



(3.2.18) 



In conclusion, [X [s) ,Y {s))t<s<T solves the following deterministic problem of ordinary (coupled) differential 
equations, almost surely. 



X(s) = f(s,X{3),Y{s)), 

Y{3) = -g{s,X{s),Y{3),Q), t<s<T, 

X{t) = x,Y{T) = h{X{T)). 

Given t,t' e [0,T], x,x' e M", consider (X|'*'^, F/'*'^, Zf'*'^)t<s<T the unique solution in Ut 



(3.2.19) 



Xe,t.x (5) ^ a; + j-^ f(r, X^^*'^ (r) , F^'*'^ (r))dr + ^e // cr(r, X^'*'^ (r) , F*^^ (r))dVK(r), 
y^'*'=" (s) == /i(X^'*'^ (T)) + /J g(r, X^^*'=" (r) , r^'*'"= (r) , Z^^*'=" (r))dr - /J Z*^'*^^ (r) dVF(r), 
t<s<T, 

extended to the whole interval [0,T], putting 

V < s < t X-^'*'^ (s) = x; y^'*'^ (s) = y^'*^^ (i) ; Z^'*'^ (s) = 0, 



(3.2.20) 



(3.2.21) 
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and consider {X^'* '^ (s) ,Y^'^ '^ (s) , Z*^'* '^ {s))t'<s<T the unique solution in A/j' 

' j^e.t'.o;' ^g) = x' + /^^ /(r, X^'*'-^' (r) , y^.*'.^' (r))dr + ^e /^^ cr(r, X^^*'^^' (r) , r^.«'.^' (r))dT4^(r), 
Ye.t',x' (-^^) ^ /^ I'^^e^t'.a;' (-y)A _^ jj ^(^^ x^^*'^^' (r) , F-^'''.^^' (r) , Z"'*''""' (r))dr - /J Z^^*''="' (r) dVK (r) , 
t' < s < T. 



extended to the whole interval [0,T], putting 



VO < s < i' X^-^'-''' is) = X- y^'*''^' is) = y^'*''^' (t') ; Z"'*''^' (s) = 



Using the estimate (3.2.9), such as in the corollary I.4 of [14], we are lead to 



E 



''^Po<s<T 



X'^*^^ (s) - X 



e.t .X 



(^) 



+ /o" 



(s)-Z'^'*'^ (s) 



supo<s<T 



„'|2 



yet,a; (-^^ _ ye 



(^^ 



<a|a;-a;'r + /3(l + |a;r)|i-i'| 



/|2 



(3.2.22) 
(3.2.23) 

(3.2.24) 



where a, /3 > are constants only depending on Ci, A. 
Finally, 



|u^ (t, a;) - u" (t', a;')r < a |a; - a;f + /9(1 + \x\^) \t - t'\" 



(3.2.25) 



which proves that (u^)e>o is a family of equicontinuous maps on every compact set of [0, T] x M". We can apply 
Arzela's Ascoli Theorem and conclude that the convergence of w^ to u, where u{t, x) = Y^'^ , is uniform in [0, T] x K, 
for every K compact subset of M'^. 

Taking the limit in e ^ in (3.2.25) we get 



\u (t', a;') - u (t, x)\ <a\x- x'\ + ^{1 + \x\^) \t - t'\" 



(3.2.26) 



for all (t,x), {t' ,x') G [0,T] x M", which proves that the function u, limit in e — > 0, is Lipschitz continuous in x 
and uniformly continuous in t. The boundedness of u is given by (2.7), since this bound is uniform in e. Using 
theorem 5.1 of [28] we deduce that u^ is a viscosity solution in [0,r] x W^ of (2.2). 

Since the coefficients of the quasilinear parabolic system are Lipschitz continuous, by the property of the compact 
uniform convergence of viscosity solutions for quasilinear parabolic equations (see ]9] for details) , we conclude that 
■u is a viscosity solution of (3.2.2). 
Moreover, let w : [0, T] x M" ^ "" ^- ~ ^^'^ 
in t for (3.2.2). Fixing {t,x) G [0,T] x M", we take the following function 



be a C;,' ([0, T], R") solution, continuous Lipschitz in x and uniformly continuous 



V- : [t,r]^R", 
V'(s) : ^v{s,X'^^{s)). 
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Computing its time derivative: 



d 



i^,, -_ |,,.....„,,^|i,,....,„2£;i(£)) 






{s, x'-^ {s)) + J2— is, x*^- is))f{s, x'-^ (s) , y*^- (.)) 



ds ' •'^-^ 9a;- 

= -g{s,X''^s),v{s,X'^^s)),0). 

As a consequence, i;(i, x) == u(t, X*'^ (i)) = u(t, x), under the hypothesis of uniqueness of solution for the system 
of ordinary differential equations (3.2.1). So, under these hypothesis, we have a uniqueness property of solution for 
(3.2.2) in the class of C\' ([0,r] x K"), which are Lipschitz continuous in x and uniformly continuous in t. D 

Remark 17. The result presented above (theorem 16) is stated under the assumptions (A. 3) and (A. 4), which 
ensure existence and uniqueness of solution of (3.2.3) for a local time T < C , where C is a certain constant only 
depending on the Lipschitz constant Ci . Under these assumptions, our results and the existence and uniqueness 
of solution for the FBSDEs (fl4j, [^3]) do not depend on results for PDE's but only on probabilistic arguments. 
However, it is possible to extend theorem 16 to global time. In order to do it and to prove the important properties 
in fl4j, results of deterministic quasilinear parabolic partial differential equations [20] are used. If we maintain the 
assumption of smoothness on the coefficients of the FBSDEs (A. 3) and (A. 4), the Four Step Scheme Methodology 
ensures existence and uniqueness of solution for (1.1), and if we remove this requirement of smoothness on the 
coefficients of the system (1.1), a regularization argument in [14] is used in order to prove it. In this setting, we 
can also conclude the claims (2) and (3) of theorem 16. 

4 Large deviation principle 

In this section, we study the Freidlin-Wentzell's large deviation principle for the laws of the family of processes 

{ (X-'*'- (•) , y-'*'- (•)) },g(o_ij e S^ (t, T- K") X S^ [t, T: M") 

as £ ^- 0. 

To begin, let us recall the following definitions mainly from [12]. 

Definition 18. If E is a complete separable metric space, then a function T defined on E is called a rate function 
if it has the following properties: 

(a) Z : _E — > [0, +oo] , X is lower semicontinuous. . 

(b) // < a < oo, then Cx (a) = {x E E : I (x) < a} is compact. 

Definition 19 (Large deviation principle). If E is a complete separable metric space, B is the Borel a-field 
on E, {/ig : e > 0} is a family of probability measures on (E, B) , and T is a function defined on E and satisfying 
(4-1), then we say that {/^el^^o satisfies a large deviation principle with rate I if: 

(a) For every open subset A of E , 
lim inf e^o e log tie [A) > - inf Z (.g) . 

(b) For every closed subset A of E, 
lim supg^o e log Me (^) < - iiif ^ (g) , 

geA 

20 



here the infimum over the empty set is defined to be +00. 

For our purpose recall the following results from [1, 3, 4] and suppose that / and a are independent of t. 
Proposition 20. Consider 

\ (7° : R" -^ R"^^''. *- ^ 

Assume that f^, a'^ satisfy Lipschitz continuous functions, with suhlinear growth. 

Consider 

/^ : M" ^ M", 

Assume that f^ , a^ are uniformly Lipschitz continuous, have suhlinear growth and converge uniformly to f := f^, 
a := a^ , respectively. Then, the family {X^'*'^ (•) : < e < 1} of random variables of the solutions of the following 
perturbed stochastic differential equations 



dX^'*'^(r) = f^X''-*'^(r))dr + y^cr^{X^^^''^{r))dW{r), 
X^'*'^ (t) = X, 



(4.5) 



obey a large deviations principle in C {[t,T] ;R"), the space of the continuous functions f : [t,T] — > R", with the 
good rate function I defined as 

I {^) = inf I ^ ll^llL(t) ■■Ven.it), ^' (s) = / {^ (,s)) + a{<l, (s)) ^' (s) 

Moreover, the level sets of T are compact and for every Borel subset A of C {[t,T] ;M") , we have 

- inf Z(5) < liminfelogP(X^'*'"= G A) < limsupelogP (X^'*'"= G A) < - inf I(.g). 

gGi s^O e^O geA 

The proof can be seen in [4, 5]. Next we introduce a very important result in Large deviation theory, to transfer 
a LDP from one topology space to another one. 

Lemma 21 (Contraction Principle Theorem 4.2.23 Page 133 in [12]). Let {lit] be a family of probability 
measures that satisfies the Large Deviation Principle with a good rate function L on a Hausdorff topological space 
S , and for e G (0, 1] , let f^.'S -^ Q be continuous functions, with {Q,d) a metric space. Assume that there exists a 
measurable map f : S ^ Q such that for every a < +cx), 

limsup sup d {fe (x) , / (x)) = 0. 

e-!-0 {x:I(x)<a} 

Then the family of probability measures {/Xe o f~^\ satisfy the LDP in Q with the good rate function 

I {y)=mi{I{x):y^f{x)}. 

For our purposes, we give the following 

Definition 22. Let t G [0, T] . We define the mapping F" ■.C{[t,T]; R") ->C{[t,T]; R") by 

F" (tp) = [s -^ u^ (s, V- (s))] , < t < s < T, V G C ([i, T] ; R") , (4.6) 

where u^ is given by (2.1). 
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As we have seen before 

From now on, for e ~ 0, u and F stand for m° and -F°. 

In order to prove the uniform convergence of the mapping F^ , we need the following formula 

||F^(^)-i^(^)||= sup \u'{s,ip{s))-uis,ipis))\, ^eC([i,r];R") 
tse[t,T] 

or 

||F^(^)-F(^)||= sup y^'«'^(^)(s)-r'^(^'(s) 

se[t,T] 

We have the following: 

Theorem 23. Assume (Al) and (A2) hold. Then, 

i) The family (X^'*'^ ('))ee(o il satisfy, as e goes to 0, a large deviation principle with a rate function 

{ve-Hl(t): 0'(s) = /(0(s),«(s,0(s)))+<T(0(s),u(s,0(s)))ip'(s), se[t,T]} ^ 

for g e C {[t,T];W-) . 

a) The family (y^^*-^ ('))ee('o il satisfy, as e goes to 0, a large deviation principle with a rate function 

l2 (0) = inf {Il (^) : F (^) (s) = (s) = u (s, ^ (s)) , s G [i, T]} . 

Proof. Indeed, noting (Al) and Lemma 14, the first assertion has been obtained in Proposition 20. We are going 
to prove the second assertion. By virtue of the contraction principle Lemma 21, we just need to show that F^ , 
e € (0, 1] are continuous and {F^} converges uniformly to F on every compact of C ([0, T] ; K") , as e tends to zero. 
Proving the continuity of F^ : 

Let £ > and x £ C([t,T],M.'^). Let (a;„)„gN be a sequence in C{[t,T],M.'^) converging to x in the uniform 
norm. Fix S > 0. Since ||a;„ — x\\^ -^ 0, there exists M > such that \\x\\^ , Ha^nlloo — ^■ 

Due to Lemma 5 and Lemma 7, we know that u^ is a continuous function in [0, T] x W^ and u'^ is uniformly 
continuous in [t,T]x K where K = B{0,M) C M''. 

There exists 77 > such that for s,si G [t,T] and z, zi G iiT | s — si |< 77 and | z — zi |< 77, we have 
\u'^{s, z) — u^(si, Zl)| < S. 

Since x„ — > a; in C([t,T],K'^), fix tiq G N such that for all n > no we have \\xn — x\\^ < rj. 

For all r G [t, T] and for all n > uq a;„(r), x(r) G K and 

|M^(r,a;(r)) - u^(r, a;„(r))| < (5. 

So, we conclude that _F^(a:„) ^- F'^{x) , which proves the continuity of F^ in the point x G C([t, T],R''). Next let 
us show the uniform convergence of the mapping F"^. Consider a compact set /C of C ([0, T] ; M") and let 

c = {^{s), (^G/c, 5G[t,r]}. 

Clearly, £ is a compact set of M". By Lemma 6, there exists a positive constant C3 such that 

sup ||F^ (ifi) -F{ip)\f = sup sup y^'''.'^('') (s) - J^^''^^") (s) ^ 

(^e/C ¥'6'CsG[t,T] 

< sup sup \Y'''-^^{s)-y'-^^is)f <C3y^. 

xeCse[t,T] 
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The proof is complete. 

Remark 24. For fully coupled FBSDEs, that is, 



D 



X^'*^^ (s) = x + // / {X^-*^^ (r) , y^'*^^ (r) , Z^-*^"" (r)) dr 

+Ve // a (X^'*'"" (r) , Y^^^'"' (r) , Z^'*'"" (r)) dW (r) , 



Y 



e,t.x 



- /J Z'^'*^^ (r) dVK (r) , < t < s < T, 



first of all, note that trajectories of the process (^Z^'*'^ {^))(t<s<T) '^'"'^ ^"^ continuous in general. As a matter of 
fact, under the assumptions of Theorem 2.6 in [30], we know only that Z^-^'''^ belongs to the space M.^ (O, T; M™^''), 
which allow us to arbitrarily change the values of the process Z^'^'^ in any P-null set. In particular, {t} x fl is a 
P-null set, which means that Z^^*'^ can be any m x d matrix. Hence, we can not get Lipschitz property for v^ . This 
issue will be carried out in our future publications. 

Remark 25. Under the conditions (A3) and (A^), we have the same conclusion as in Theorem 23. The proof is 
similar we omit it. 



5 Appendix 

Proof of Lemma 5. 

Proof. Applying Ito's formula to \Y^'*'^ (•) — Y^''^'^ {■)\ on [s,T] , we have 



lY'^^''' (s) - Y^'^^y {s)f + 



T-J's 



(\Z'-*^'' (r) - Z'-'^y {r)\^) dr 



= E- 



^a 



\h{X^'*^'={T)) - h {x^'^^y {T))\ 



-2E 



J^s 



(ye,t,x (^^ _ ye,t,y ^^^^ (^je.t,x ^^^ _ fe,Uy (^)^ j^ 



(5.1) 
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which yields 



E 



sup lY^'^'^s) - Y'^'-y {s)\ 

t<s<T 



Ie 

2 



(|Z"'*'^(r)-Z"'*'2'(r)|^)dr 



< 



Ci jX-^'*'^ (r) - X-^'*'^ (r) 1^ + (3Ci + 2Ci2) JF^'*'^ (r) - F^'*''' (r) |^ dr 



+E 



< CfE |X='*'"=(r)-X^'*'^(T)| 



+E 



(Ci |X^^*'^ (r) - X-^'*^^ (r)l^ + (3Ci + 2Ci2) jl"^'*^^ (r) - r^'*'« (r)|^) dr 



< (3Ci+2Ci)|e |X^'*'"=(r)-X'^'*^«(T)| 



+E 



f ("jX"'*'^ (r) - X'^'*^?' (r)|^ + jy^'*^^ (r) - y^'*^^ (r)f ) dr 



Applying Ito's formula to X*^'*'^ (s) — X*^'*'** (s) yields that 

= I 2(X"'*'^(r)-X"'*'2'(r)) (.f'*^^(r)-,f'*'2'(r))dr + e /" |cr"'*'^ (r) - fT^'*'« (r)|^ dr 
+\/i / 2 (X^^*'^ (r) - X-^'*^^ (r)) {<t'^''^ (r) - ct^^*'^ (r)) dW (r) . 



By Burkholder-Davis-Gundy's inequality, there is a constant C3 > such that 

E sup |X"'*'"(s)-X^^*'^(s)| 

t<s<T 



dr 



< C3<;(l + Vi)E / (|X^^*'"(r)-X^^*'^(r)|')dr +E / ([.f^*'" (r) - /"'*^M0|') dr 

((eCi^ + 2) jX-^'*^^ (r) - X^^*'2' (r)f + QCl |y"^*'=" (r) - F^^*'^ (r)^ dr 
(6(72 + 2) / f|X'^'*^^(r)-X^'*'^(r)hdr + 6rC2 sup IF-^'*'^ (s) - F^^*'^ (s)| 



< C3E 



< C3E 



(5.2) 



(5.3) 



(5.4) 
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Applying Ito's formula to (X*^'*'^ (r) — X^-'^^y (r)) (y^^*'^ (r) — Y^'''-'^ (r)), and by virtue of assumption (A2), we have 

C2E Hx"'*'^ (T) - x^-^^y (T)f 1 - E [(x - y) (y^'*'^ [t) - Y'^^^-y {t))\ 

< E [/i (X-^'*'^ (T)) - /i (X^'*'^ (T)) (X-^'*'^ (T) - X^'*'^ (T))] 

(|X^^*'=" (r) - X^'*'^ (r)l^ + |y^'*'^ (r) - y^'*'^ (r)|^) dr 



< -C2 (1 + V^) E 



< -CoE 



The inequality (5.5) can be rewritten as follows 

C2E[|A^'*'"=(r)-A:'^'*^^(T)| 



-C2E / |X'^'*^^(r)-X"'*'2'(r)|V|y"'*'^(r)-y"^*'^(r)| 



dr 



< 



[(.T - y) (y-'*-- (t) - y-'*'« (t))] . 



Combining (5.5) and (5.6), we have 



sup |y"'*'^(s)-y^'*'«(s)| 

t<s<T 



< 



E 

(3Ci + 2Cf 



-E 



\Z'-t^^ (r) - Z'-'^y {r)\ dr 



Co 



[{x v) (y-^*'- (t) - y^^*'^ (i))] 



< i|yM,-(i)-yM^2'(i)f + l('-^^i + 2C2^ , ,2 



C2 



\x-v\ 



Therefore, 



E 



sup |y^'*'^(s)-y"'*'2'(s)| 

t<s<T 



•3Ci + 2C?\ 2 

' — c^ — ) l^^^l 



By Gronwall's inequality, there exists a positive constant C4, such that 



E 



sup \X^''^'' {s) - x^'^^y {s)\ 

t<s<T 



<Ci\x-y\ , 



(5.5) 



(5.6) 



(5.7) 
(5.8) 

(5.9) 



where C4 depends on Ci , C2, C3, and T. Finally, taking Ci — max < I — ^-^ — ~ ) 1 ^^4 f" 1 we get the desired result. D 
Proof of Lemma 6 
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Proof. Applying Ito's formula to |X^'*'^ (s)| yields that 






By Burkholder-Davis-Gundy's inequality, there is a constant C5 > such that 



E 
< C5 
= Cs 



sup \X'-''^{s) 

t<s<T 



(|X^'*'^(r)hdr+ /" (|/'^'*'^(r)f)dr + eE f (icr-^'*'"^ (r)^ 



dr 



I^M,. (^)|- ^ ifs.t,. (^) _ ^^ (^^0^0) ^ /(r,0,0)| 



+e jcr^'*'"' (r) - a (r, 0, 0) + cr (r, 0, 0)|^ dr | 



< C5 



(|X^'*'" (r)l' + eC? (|X^'*'- (r)l' + |y^'*'" (r)l') + 2 |/ (r,0,0)|') 



/7 
(|f7(r,0,0)hdr 



dr 



+4eCfl 

< C5{2 + 8Cf) 

< C5 (2 + SC'l) E 
+ / (|/(r,0,0)|V|<T(r,0,0)|')dr 



(|X^'*'" (r) f + |r^'*'" (r)l' + 1/ (r, 0, 0)|' + \a (r, 0, 0)|') dr 



X"'*'^(r)rdr + T sup |r"'*'^(s)| 

t<s<T 



Second applying Ito's formula to |y^'*^^ (•)! on [s,T] , we have 



|y"^*'="(s)f + E 



^a 



(|Z"'*'^(r)f)dr 



7J^s 



|/i(x'^'*'^(r))| 



r^. 



ye,t,x (^) je^x ^^^ j^ 



(5.10) 



(5.11) 



(5.12) 
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which yields 





E 


< 


E 


< 


E 



sup \Y''*-^s) 

t<s<T 



Ie 

2 






dr 



h{X'^*'^T))-h{0) + h{0)\ 



(|r"'*'^ (r)|^ + 1/='*'^ (r) - / (r, 0, 0, 0) + / (r, 0, 0, 0)|^) 



dr 



< E 



2Ct\X''-'''''{T)\ +2\h{0)y 






< 



/((^+3C?)|y-(r 

+ r(|/(r,0,0,0)p)dr+|/i(0)|4. 



dr 



Set 



Q^(3C? + 3+^)(|M0)|V/' 



|g(r,0,0,0)rjdrl. 

Applying Ito's formula to X^'*'^ (r) y^^*'=^ (r), and by virtue of assumption (A2), we have 
E [X^^*'"" (T) /i"^*'^ (X^'*'^ (T))] - E [xY^^*'"" (t)] 

< -C2(l + Vi) 

+E 

< -^E 

2 



(|X"'*'^(r)|V|r"^*'="(r)|^)dr 
(X^^*'" (r) .9 (r, 0, 0, 0) + F^'*'" (r) / (r, 0, 0) + ^iZ^^*'" (r) ct (r, 0, 0)) dr 
(|X"^*'="(r)f + |y^'*'^(r)|^)dr 



-E 



1 

^2C2 

< -^E 
2 



|5(r,0,0,0)|' + |/(r,0,0)ndr 



ViZ^^*'="(r)a(r,0,0)dr 



hx^^*'"" {r)f + \Y^'*^'= {r)f\ dr 



Ml, 



where 



Ml = 



/; |.g(r, 0,0, 0)1 V|/(r, 0,0)1' dr 



2C2 



(^/iZ'^■*^^(r)a(r,0,0))dr 



(5.13) 



(5.14) 
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On the other hand. 



> E 

> E 



[X^'*^"" {T) h^'*^'= {X^'*''-" (T))] 

^(^X^Ax (y) _ 0) {h^^*''' (X'^'*'^ (T)) - h^^*''-" (0) + /i'^'*^^ (0))] 

C2 IX^^*'"" {T)f + X^^*''-" (T) h^'*^'= (0) 
C2 IX^'*'"' (T) 



e.M,^M2 |X-^*'^(r)|' l/i- 



(0)r 



2a 



where a > large enough such that C2 — ^ > 0. 
Then, we have 



Co 



E 






< E [xY"-*''' (t)] + Ml + a 



|/je,t.^(0)|^ 



Setting 



we have 



|/i^'''="(0)|'' 



M2 = Mi+a 2 



E 






< 



E [a;y^^*'^ (t)] + M2 



Noting (5.13), we obtain that 



E 



sup |r^'*'''(s)P + i / |Z'^'*'^(r)|^dr 

t<s<T 2 7j 



< (3C?+3+^ 



[.xy^'*'^ (i)] + M2 
c 



Ce 



< 



|yM,.(^)|2 \:,\^{iCl+i + ^y 



2 2C2 

SCf + 3 + -^) M2 



C 



+ C6. 



(5.15) 



(5.16) 



(5.17) 



(5.18) 
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Define 



The expression (5.18) can be rewritten as 



M 



iCi + 3 + g^ 



6Cf 



C 



E 



sup |r^>*^^(s)P + i / |Z"^*'="(r)|^dr 

t<s<T 2 

,t,x /^\|2 I |2 ,->2 



< 



|y^'*'^(i)r \x\^M 



+M a- 



\h^''--{0)f /o (|5 (r, 0, 0, 0)1^ + |/(r, 0,0)1' 



dr 



2^2 



+M2 / |cr(r,0,0)rdr + 



Z^'''^(r) dr 



Cfi 



Consequently, we get 



(5.19) 



E 



sup |y"'*'^(s)| 

t<s<T 



|Z"'*'^(r)f dr 



2 ^2 ~ / |/i^'*'nO)|' /o ( l5(r, 0,0, 0)1^ + 1/ (r, 0,0)1' 



dr 



2^2 



+2M-' / |cr(r,0,0)rdr + 2C6 



< max 



{E,Af^}(l + N^) 



where 



(5.20) 



E = M a- 



/iM,.(o)|2 /; (Iff (r, 0, 0, 0)1^ + |/(r, 0,0)1' 



dr 



2^2 



+ 2M |a(r,0,0)rdr | +2Ce 



By Gronwall's inequality, we derive 



E 



sup |X^'*'^(s)| 

t<s<T 



<C7 i+ixn. 



(5.21) 



where C7 is independent of e. Taking 



Co = max 



|max|s,M2|,C7} 



we get the desired result. 
The proof of Lemma 7. 



D 
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Proof. Suppose that ti > ^2- By standard estimates and Ito isometry we have 

E 



sup |X^'*i^^(s)-X"'*^^^(s)| 

.ti<s<T 



< 2{ti-t2)l 

< iti-t2)E 

r-T 



/(r,X^^*^'^(r),y^'*^'^(r))| dr 



cr(r,X'^'*^'^(r),y^'*^'^(r))| dr 



16TC2 f sup |X'^'*^^^(r)|V sup JF^^*^^^ (r)|^) + 2 / |/ (r, 0, 0)|^ dr 

\t2<r<ti t2<r<ti ) Jo 



+2 / |cr(r,0,0)rdr 
Jo 

It follows from Lemma 6 that 



E 



sup |X^'*i'^(s)-X^'*^'^(s)| 

ti<s<T 



< (ti-i2)e, 



where 



e 



lerC^Ca (l + |a;|') + 2 / |/ (r, 0,0)|' dr + 2 / |a (r, 0,0)|' dr 



Similarly as in Lemma 5, we have 



E 



sup |y"'*i'^(s)-r^'*2'=^(s)| 

tx<s<T 



|Z^'*i'="(r) -Z^^*^'"=(r)r dr 



< 



(3Ci + 2Cl) |e |X'^>*i>=^ (T) - X^^*^'"= (r)| 



+E / |X"^*i'^ (r) - JC^'*^'^ (r)| + |y^'*i>=^ (r) - F"^*^'^ (r)| dr ^ . 



By Gronwall's inequality, we obtain 



and 



E 



sup |y^'*i'"(s)-y^^*^'"(s)|' 

Lti<s<T 

< (ti-i2)e^(3Ci + 2Ci2)(l + T)e 



Z'^''i'^(r)-Z'^''^'^(r) dr 



Set 



and 



< 2 (ti - t2) (3Ci + 2C;i2) e (1 + T + e^ (3Ci + 2Cl) (1 + T)) 



Cg = e^ (3Ci + 2Cl) (1 + T) e 

Ce = 2 (3Ci + 2C2) 6 (l + T + e^ (3Ci + 2^2) (1 + T)) 

C3 = max{e,C5,C6}. 



We get the desired result. 



D 
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The proof of Lemma 8. 

Proof. Analogously, applying Ito's formula to jy^i'*'^ (.) — ye2,t,2: ^.-jj q^ [t,T] , by the method used above, we 
have 



E 
< E 



sup |r^i'*'^(s)-r^^'*'^(s)|^ 

t<s<T 



t,x /„\|2 



Z'^-'''^ {r) - Z"^-'''' {r)\ dr 



(C'l |x^i'*'^ (r) - X-^^'*'^ (r)l^ + (3Ci + 2Cf) \y^-^-*^^ (r) - y^^'*'^ (r)^ dr 



< Ci^E |X'^i'*'^(T)-X'^2'*'^(r)| 



< 



(Ci |X"i'*'^ (r) - X-^^'*'^ (r)|^ + (3Ci + 2C^) Ir-^i'*-^ (r) - r-^^-*'^ (r)|^) dr 
(3Ci + 2Cf) |e |x^1'*'^ (T) - X-^^.t.o: (-j,^|2 
+E /" f jX-^i'*'^ (r) - X^='*'^ (r)|^ + |r=i'*'^ (r) - F"^'*'^ (r)f ) dr i . 



(5.22) 



Applying Ito's formula to X^^'^'^ (s) — X^^'*'^ (s) yields that 

|a:"1'*'="(s)-a:"^^*'="(s)|^ 

= / 2(A:"i^*'="(r)-X'^^'*^^(r)) (/"i^*'="(r)-/^^^*'="(r))dr + e / |ct^i^*'=" (r) - tr'^^'*'^ (r)|^ dr 
+Vi / 2 (X^i^*>=^ (r) - X-^^'*^^ (r)) (cr^i^*'=" (r) - a^^^*'=" (r)) dW^ (r) . 



By Burkholder-Davis-Gundy's inequality, there is a constant Cg > such that 



E 
< Cs 



sup X^i'^'" (s) - X^^^''" (s) 

t<s<T 



(l + \/i) |X"i^*'="(r)-X'^^'*^^(r)f dr 



-E 



I -eEi,t,x f^\ __ f£2.*-^ (j-)! (Jt- 



(Vi + e)E 



|cr"i'*'^(r)-cr"^'*^^(r)f dr 



< CsE 



< CsE 



((eCi^ + 2) jX-^i'*^^ (r) - X^^'*'^ (r)|^ + eCi^ |r^i'*'^ (r) - y^^'*'^ (r)l^) dr 



(6Cf + 2) |X'^i'*'^(r)-X"^'*'^(r)fdr + 6rCi2 sup JF^i'*'^ (s) - F"^'*'^ (s)| 

Jt t<s<T 



(5.23) 



(5.24) 
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Once again applying Ito's formula to (^X^^'^'^ (r) — X^^'*'^ (r)) (y^i^*'^ (r) — y^a.t.a; ^^^-j_ ,^^^ by virtue of assumption 
(A2), we have 



C. 



-Y'ei,t,x /rp\ -Ye2,t,x /rp\ 



< E [/i [X^^^^'^ (T)) - h [X"^'*^'-' (T)) (X'^i>*^^ (T) - X^^^''"" (T))] 



< -CoE 



+ 2(^/^-V^)E 



< -C, 



f f |X"i^*'^ (r) - JC^^'*^^ (r)l^ + |y^i^*.=^ (r) - y^^'*'^ (r)|^) dr 
cr"^'*'^ (r) (Z"i^*'=" (r) - Z"^'*'^ (r)) dr 

(\X'^ut.x (^^ _ ^e2,t,x (^^|2 _^ |^ei.t,a= ^^^ _ ye2,t,x (^^|2\ ^^ 

(|(7"^'*'" (r) -CT(r,0,0) +a(r,0,0)|' + |Z"i'*'" (r) - Z"^'*'" (r)|') dr 
(|X"i'*'^ (r) - X-^^'*'^ (r)l^ + jy^i'*'^ (r) - F"^'*'^ (r)|^) dr 



+ (a/si - \/e2) 



4TCM sup |X'^^'*^"'(r)f + sup ly^^.t.a: (-^^ 

t<r<T t<r<T 



(|Z"i'*'^(r)|V|Z'^^'*'^(r)|^)dr 



+ 2^ |cr(r,0,0)rdr 



(5.25) 



From Lemma 6, we know that 



( sup |X"^'*'^(r)P+ sup |y"^'*'^(r)P+ / IZ^^'*'^ (r)P (r) + / IZ^^'*'^ (r)P (r) dr 

\ t<r<T t<r<T Jt Jt J 



< 4C2[l+\x\ 
Hence, combining (5.22), (5.25), (Al) and (A2), we have 



E 



sup |y'^i'*'^(s)-y'^2,t,:r(^^| 

t<s<T 



1. 



|Z'^i'*'^(r)-Z^^^*'^(r)rdr 



< ^^'-^K ,{3C,+2C!), 



(5.26) 



where 



M. 



(3Ci + 2Ci2) [(l + |a;|') (2C2 + STCfCa) + 2 /g^ |f7 (r, 0, 0)|' dr 



Co 



By Gronwall's inequality, we also obtain 



sup IX"'-'^'' (s) ^ X"^^'-'' (s) 

t<s<T 



<C9(V^-V^) 
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where Cg is independent of e. Taking 



f^ M3(3Ci+2C2) 

1 ^9' 7^ 



d = maxs „„, 



we complete the proof. D 
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